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A generalized fluctuation-response relation is found for thermal systems driven out of equilibrium. 
Its derivation is independent of many details of the dynamics, which is only required to be first- 
order. The result gives a correction to the equilibrium fluctuation-dissipation theorem, in terms of 
the correlation between observable and excess in dynamical activity caused by the perturbation. 
Previous approaches to this problem are recovered and extended in a unifying scheme. 
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The fluctuation-dissipation theorem is a standard 
chapter in statistical mechanics P, S, A system m 
thermal equilibrium has statistical fluctuations propor- 
tional to its response to external perturbations: a small 
impulse changing the potential U U — hgV at time 
s, will produce a response s) = S{Q{t)) /Shs in a 

quantity Q at time i > s of the form 



R''^vit^s)=P-{V{s)Q{t)) 



cq 



(1) 



where {V{s) Q{t))^^ is a correlation function quantifying 
the equilibrium fluctuations in absence of any pertur- 
bation, and the proportionality constant /3 = l/ksT is 
the inverse temperature. An early example of this theo- 
rem is present in Einstein's treatment of Brownian mo- 
tion, where the diffusion constant, expressed as a veloc- 
ity autocorrelation function, is found proportional to the 
mobility. Other famous examples include the Johnson- 
Nyquist formula for electronic white noise and the On- 
sager reciprocity for linear response coefficients. 

So far, approaches deriving a fluctuation-dissipation 
relation (FDR) for nonequilibrium 0, d, 0, 0, i, E E, 
111 . Il2l . [l3| have not found a physical unification and do 
not appear as textbook material. One reason may be that 
previous work has not been seen to identify a sufficiently 
general structure with a clear corresponding statistical 
thermodynamic interpretation. Today, such an interpre- 
tation has become available from advances in dynamical 
fluctuation theory for nonequilibrium systems. 

Aiming to provide a simple and general approach to 
FDRs, in this Letter we put forward a FDR for nonequi- 
librium regimes in a framework that may represent a uni- 
fying scheme for previous formulations. Our main result 
can be found in a general formula, Eq. ([6]) below, which 
can also sometimes be rewritten as ([7]) or ([TT|). 

In order to go beyond equilibrium and beyond formal 
perturbation theory, it is important to recognize in the 
right-hand side of ^ the role of entropy production, 
as usual governing close-to-equilibrium considerations. 
Eq. ([1]) expresses the correlation between the dissipa- 
tion represented by entropy production (energy change 
divided by temperature) and the observable Q{t). Here, 
for perturbations of a nonequilibrium system (that has 



already a non- vanishing set of flows and hence a non-zero 
entropy production) we will rather speak of the excess of 
entropy produced by the perturbation hs V. 

The second ingredient that is essential in nonequilib- 
rium is still a less known quantity, called traffic jl4| or 
dynamical activity [l5l |. first introduced in [iq| . Perhaps 
activity has been somewhat overlooked in the past be- 
cause it plays a truly significant role only beyond linear 
order around equilibrium [ij, [l7| , and it is not part of the 
picture in standard irreversible thermodynamics. Our 
present results will in fact further clarify the relevance of 
dynamical activity. To get a first idea of its meaning, we 
note that activity often measures the frequency of tran- 
sitions in a trajectory: this is a time-symmetric aspect of 
dynamical fluctuations [lij because it does not depend on 
the direction chosen to span the trajectory. On the other 
hand, it is well known that entropy production is time- 
antisymmetric [lit I changing sign upon reversal of time 
(fluxes are inverted if time is run backward). Thus, both 
quantities arise naturally as two complementary pieces of 
the space-time action, as reminded below. 

We consider a system in an environment at inverse 
temperature (3 (with ks = 1), where a time-independent 
nonequilibrium condition can be imposed by external 
driving fields or by installing mechanical displacements 
or chemical gradients at the boundaries of the sys- 
tem. Moreover, the system evolves according to a 
Markovian/first-order dynamics: for any two observ- 
ables / and their correlations at times s < t satisfy 
^{f{xt)g{xs)) = {{Lf){xt) g{xs)) where L is a linear 
operator, called generator, acting on observables, and xt 
is the state at time t at some reduced, e.g. mesoscopic 
level of description: e*^/(a;) = {f{xt)\xo = x), t >Q (in 
the sequel we abbreviate f(xt) = f{t)). The generator 
i is a standard tool: for jump processes it is the ma- 
trix with transition rates as off-diagonal elements, while 
L = force • V -I- A//3 for overdamped diffusions. How- 
ever, it is our aim to leave the generator L as unspecified 
as possible, to emphasize the generality of the results, 
referring e.g. to Q for further background. 

Our formulation is based on the point of view that 
the system can be described by a distribution 'P{u3) over 
its possible trajectories u) (space-time paths), giving the 



2 



state LOs of the system at each time s < t. In general 
'P{<jl3) = ^(wlwo) /^("^o); namely ■p(ci;) is the probabil- 
ity of the initial state ujq, i.e. ^(wq), times the probabil- 
ity of the trajectory V(y3\ujo) given the starting point. 
The ensemble averages of a quantity 0{ui) are then 
{0{u})) = V{u})0{u>) where is a simple notation 
including also a continuum of trajectories. Here we are 
interested in comparing probabilities of trajectories for 
hs^Q (for s > 0), V^iuj) ^ P''(uj|u;o)/u(wo), with those 
of the unperturbed system. In analogy with the h = 
case, we use the notation {0{u}))^ = V^{u})0{u}) for 
averages in the perturbed system. We start by writing 

V'^iui) ^ e-^('*')p(u;) (2) 

and by focusing on the action A{Ld) ~ — \\Y\P^{ui)/V{ijJ)\. 
In this set-up we put ourselves in line with the Onsager- 
Machlup approach but outside equilibrium, possibly be- 
yond quadratic or diffusive approximations and also deal- 
ing with jump processes. 

We ask what determines the path-space distribution 
and the space-time local action governing it. To find 
an answer, it is useful to decompose Aioj) = {T{u)) ~ 
S{u}))/2 in terms of its time-antisymmetric component 
S{ijj) = A{9uj) — A{ui) and the time-symmetric T{ui) = 
A{9<jj) + A{u)) , where a time-reversed state {9u))s is equal 
to uJt-s (with reversed momenta when applicable). In the 
time interval [0,t], in standard physical situations one 
finds that 



htV{t) - hoV{Q) 



d.s— — V(s) 
ds 



(3) 



is the excess entropy flux from the system to its envi- 
ronment. Excess is always meant in the sense of the 
perturbed process with respect to the original one. From 
basic thermodynamics we know that this entropy flux is 
/3 Q, where Q is the heat flowing to the thermal environ- 
ment, that is: minus the change in potential energy of the 
system minus the work performed on the system. The 
reason for ^ is the physical condition of local detailed 
balance, which ensures that the ratio of probabilities of 
the forward with respect to the backward trajectory is 
given by the exponential of the entropy flux in the for- 
ward trajectory [isj . 

The time-symmetric part T{lS) of the action A is the 
excess in activity, having an essential role in dynamical 
fluctuation theory for the large deviations of the occu- 
patioiis, and in nonequilibrium studies of phase transi- 
tions [ij, [l^. Being in linear response theory (small /is), 
we are actually more concerned with the first order in 
excess activity r(a;, s) = j|-T(u;) j/j^^o- The excess in 
activity quantifies how "frenetic" is the motion in the 
perturbed process, compared with the unperturbed one. 
For example, for a Markov jump process with transition 
rates W[x y) between states x ^ y, the traffic is equal 
to twice the time-integrated escape rates over a trajec- 
tory, 2 ds '^{'^s y) and its excess r(s) (here 



t(ci;, s) = t(s) depends only on the state uJs) to linear 
order in V is 



r(s) 



2^M/(, 

y 



^§lViy)-V{u.,)] 



y 



-y){ 

Wiujs ^ y)[Viy) 



1 



(4) 



dV 



The last expression is of course formal, the mathematical 
meaning being dV^/d.s = LV{s). 

Since S{u^) is already linear in hs, we have 



-A{uj)^^S{u:)-lJ^dsh, 



-{u^,s) + 0{hl) (5) 



We are now ready to derive the FDR for a single- 
time observable 0{u)) = Q{t) at time t > 0. Its ex- 
pectation in the linear response regime is {Q{t))^ ~ 
{Q{t)) + Jf^ ds hsRQvit, s), which we can rewrite with 
@ as (Q(0e-'4{-)) - {Q{t)) ^ -{Q{t)A{Lj)) = 

dshgRgvit, s). With ([3]) and ([5]), it is thus straight- 
forward to see that the response function is equal to 

RQvit, s) = f ^ {Vis)Qit)} - \ {t{u,, s)Q{t)) (6) 

which, as we will see, is a nonequilibrium FDR that gen- 
eralizes and unifies previous formulations. We stress that 
the derivation of ^ is independent of many details of 
the dynamics and hence of its generator L: we went be- 
yond formalities of a perturbative first-order calculation 
by specifying the origin of the two parts on the right- 
hand side of ^ in terms of (i) the entropy flux and (ii) 
the dynamical activity 14, 15, 1^, in excess by the time- 
dependent perturbation. 

When r(w, s) = t(s). as in ^ or for overdamped dif- 
fusions [19), formula ^ simplifies to 

RQv{t, = f ^ {V{s) Q{t)) - ^ (^^Vis) Q(t)^ (7) 

The right-hand side of ([7]) is nonzero because of the cor- 
relation of V{s) with Q{t) ioTt> s^. 

When the unperturbed dynamics is detailed balanced, 
we can recover ([T|) as a special case of ([7]) by using 
time-reversal symmetry: for s < t, (dV/ds Q{t))eq = 

where the last 



{dV/dtQis)), 



d/dt{Vit)Q{s)), 



equality is a consequence of the Markov property of the 
dynamics. Because of stationarity this is then also equal 
to — d/ds {V{s) Q{t))cq, so that in equilibrium the second 
term of ([7]) gives exactly the same contribution as its first 
term, and ([1]) is recovered. 

As an illustration of the general formula, consider a 
system with a flow of interacting particles due to the 
presence of reservoirs at different chemical potentials. We 
choose both Q ^ V = TV, where J\f is the particle num- 
ber. Hence, hg induces a global shift in chemical poten- 
tials, and the linearized excess activity r(s) = /JdjV /ds = 
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FIG. 1: (Color online) (a) Plots of x(i)//9 vs [C{t) + D{t)]/2 
for times t = 5, 10, . . . , 50 (other parameters: see legend, and 
h = -0.1, di-d„^ 0.8). (b) C{t = 50) and D{t = 50) as a 
function of the unbalance between reservoir densities, for /3 = 
1, n = 20. At equilibrium (di = d„) one has C(50) = D(50) 
as expected, while the difference between C(50) and D(50) is 
not negligible for strong nonequilibrium. 



/3i7(s) is proportional to the systematic particle current 
J{s) into the system via its boundaries, i.e. to the ex- 
pected current given the configuration at time s. Numer- 
ically it is convenient to test the time-integrated version 
of ([7]) for constant perturbation hg ~ h ior s > 0, 



x{t)/l3^[C{t)+D{t)]/2 



(8) 



with susceptibility x{t) = {{-^{t))^ ^ correlation 
function C{t) = (7V^) - {Af {0)N {t)) , and nonequilibrium 
term D[t) = - f^ds{J{s)M{t)). Again, © holds irre- 
spective of the details of interaction or driving. 

Here we simulate a boundary driven one-dimensional 
exclusion process, a simple paradigmatic model for trans- 
port. A state X is an array with empty sites a;* = and 
particles = 1. Nonequilibrium is imposed at boundary 
sites i = \ and i — n, which arc in contact with reser- 
voirs with particle densities di and d„ = 1 — di, respec- 
tively. There is a nearest neighbor attraction with en- 
ergy H{x) = — Y^"^i ic'a;'"'"^ determining the Kawasaki 
dynamics in the bulk. In this context the systematic 
current arises from two possible transitions: x y, with 

^ 1 — x^, and x with ~ 1 — x"'. Hence 

J{x) = [J\f{y)-J\f{x)]W{x^y) + [Af (z) - Af {x)]W {x ^ z) 

where e.g. W{x ^ y) oc d^e^'^I^^^^^-^^^^l/^ if a particle 
enters a.t i — 1. Examples in Fig. [Ija) visualize Eq. ([5]), 
while the magnitude of C{t = 50) and D{t = 50) is shown 
in Fig. [Hb) as a function of the driving di — d„. Far 
from equilibrium C(t) ^ D{t) and thus the correlation 
function C{t) alone cannot provide a good estimate of 
the susccptibihty x(0- 

In order to connect Eq. ([6|) [or ([7|)] with some previous 
formulations, it is easy enough to rewrite it in another 
form. First remember that causality implies RQv{s,t) = 



for s < t. Hence, for ([6]) to be valid for all s, t one needs 
/3^(F(t)Q(s)) = (T(c.,t)Q(s)) fors<t, (9) 

which is automatically satisfied in our framework: it fol- 
lows from (e^'^('^^) = 1, equivalent to {A) = to first 
order in hs, via partial integration of ([3]). By subtracting 
dHl) from (l6l) one finds 



RQv{t,s)^^(^-^^{Vis)Q{t)) 
1 



--(r(a;,s)Q(i) 



TiLu,t)Qis)) 



(10) 



In stationarity the second line in pop is an extra term 
compared to the equilibrium version ([T]). in which the 
time-antisymmetric correlation between activity r and 
the observable Q is evaluated. This term has been called 
asymmetry in studies of overdamped Langevin equa- 
tions [1] and of discrete stochastic systems 0, [IS] ■ 

Again, when r(w, s) = t(s), as in (j4|) and ([7]), equation 
^ implies that t(s) = l3LV{s). Hence ^ turns into 

/3 d P 



RQvit,s) 



{LVis)Qit)) (11) 



which is the rigorous version of ([7]). Let us now make 
contact with some previous formulations for stationary 
processes, where £ {V{s) Q{t)) = £ (1/(0) Q{t - s)) = 
— ■^{V{s)Q{t)). From our definitions, the last term 
equals —{V{s)LQ{t)). Alternatively we can think of 
Q{t) correlated with V{s) evolved backward in time by 
the adjoint L* (which generates the time-reversed pro- 
cess): {V{s)LQ{t)) = {L*V{s)Q{t)). Eq. ^ can then 
be rewritten as 

RQvit,s) = -^{L*Vis)Q{t) + LVis)Q{t)) (12) 

which is a generalization of Eq. (2.15) in Q. Indeed, in 
the context of fluctuation theory around diffusive scaling 
limits, the L* is referred to as the adjoint hydrodynamics 
in the infinite-dimensional treatment in Section 2.3 of 0]. 
It is emphasized there that the response for the adjoint 
process is typically along the reversed trajectory from a 
spontaneous fluctuation in the original dynamics. This is 
also explicit in ()12p . as the response in the time-reversed 
process amounts to the same expression upon exchanging 
L with L* and t with s. Note however that one needs 
the stationary distribution to know L*. The knowledge 
of this distribution enables also the derivation of other 
FDRs [3,i,S,[i,[2i|. On the other hand, Eq. ^ does 
not involve the stationary law, except for the statistical 
averaging. 

We also recover the interpretation of jTl|: by using 
(fT2|) . we rewrite the stationary version of ([6]) as 



i?Qy(i,s) = 



f3-{Vis)Qit)) 



f3 



{[{L-L*)V]is)Q{t)) 
(13) 
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Note that the response regains its equihbrium form [only 
the first term in the right-hand side of (|13|) ] whenever 

is a time-direction neutral observable, in the precise 
sense that LV = L*V. Furthermore one can check (say, 
for overdamped diffusions and for jump processes) that 
L — L* = where j is the probability current and p 

is the stationary distribution on the states. Their ratio is 
a drift velocity v = j / p. Hence, for a time-direction neu- 
tral potential V, the probability current is orthogonal to 
its gradient, v ■ W = 0, and the equilibrium form of the 
response is obtained, in agreement with the observations 
in [2^ . The same effect is achieved by describing the sys- 
tem in the Lagrangian frame moving with drift velocity 
V. The second term in (|13p vanishes also in this case, 
an d (11^ yields exactly the interpretation (and E q 4 5) 
of , which was recently experimentally verified [23| . 

To finish, we revisit the relation with dissipation to 
confirm the prediction in that the usual equilibrium 
relation between response and dissipation is preserved 
when taking into account only the excess heating and 
ignoring the housekeeping heat. During a specific trajec- 
tory of the system, the heat dissipation can be split into 
two parts: Qhk + Qex [2J]. The first term is the house- 



keeping heat, which is the heat produced in the unper- 
turbed dynamics. This is the heat that drives (or is the 
result of) the system out of equilibrium. The second part: 
Qexi is the extra heat generated through the perturba- 
tion, and is just 1/P time the excess entropy production, 
which we already encountered in our calculations: 



Qe 



4 (SHf 



ds 



^ {V{s)t 



ds 



(14) 



This expression is of the same form as in equilibrium, 
so we will not repeat the calculations here, only the 
conclusion: the (excess) heat is given by the imaginary 
part of the Fourier transform of Rvv, or equivalently, 
the Fourier transform of the time-antisymmetric part of 
Ryv- Moreover, we have already written down this time- 
antisymmetric part: it is the right-hand side of (jlO[) . ex- 
tended io t < s. 

In conclusion, from physical constraints on the prob- 
ability of trajectories, we have obtained a general FDR 
for the response of a driven system to the addition of a 
potential. It lifts the nonequilibrium FDR beyond for- 
mal first order perturbation theory applied to a specific 
dynamics, by identifying in general physical terms the 
statistical quantities that determine the response. Previ- 
ous formulations are recovered as specific cases. Finally, 
we have given observational significance to the notion of 
dynamical activity, which so far has mostly appeared as 



a theoretical concept in fluctuation theories. 
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